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Abstract. Diffusive shock acceleration is the prime candidate for efficient acceleration of 
cosmic rays. Galactic cosmic rays are believed to originate predominantly from this process 
in supernova remnant shock waves. Confinement of the cosmic rays in the shock region is 
key in making the mechanism effective. It has been known that on small scales (smaller 
than the typical gyroradius) high-amplitude non-resonant instabilities arise due to cosmic 
ray streaming ahead of the shock. For the efficiency of scattering of the highest energy 
cosmic rays it is of interest to determine the type of instabilities that act on longer length 
scales, i.e. larger than the cosmic ray gyroradius. We will present the results of our analysis 
of an instability that acts in this regime and will discuss its driving mechanism and typical 
growth times. 
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1. Introduction 

Galactic cosmic rays are widely believed to 
be accelerated, be it predominantly or not, in 
supernova remnants (SNRs). If solely SNRs 
are responsible for the acceleration, a signifi- 
cant fraction of kinetic energy has to be trans- 
ferred onto cosmic rays. Additionally, the ac- 
celeration needs to be efficient enough to ar- 
rive at energies typical of the knee in the spec- 
trum, being around lO'^ eV. In order to get 
those acceleration rates, Bohm diffusion is nor- 
mally assumed, since this yields the quickest 
isotropisation needed for the confinement of 
cosmic rays around the shock region. In real- 
ity, the diffusion is dependent on the power in 
magnetic fluctuations on the scale of the gy- 
roradius. Scattering of the highest-energy cos- 
mic rays requires a significant component of 
magnetic field fluctuations on the order of par- 
sec scales. A cosmic -ray streaming instability 



that initialises growth on those scales is de- 
rived from the Boltzmann equation for the cos- 
mic rays, coupled to the MHD equations of the 
background plasma. Typical growth times are 
presented. 

2. Cosmic ray streaming instability in 
the long-wavelength limit 

It has long been known that cosmic rays can act 
to trigger magnetic fluctuations. The first cal- 
culations on this were in the resonant regime, 
where energy is transferred from the cosmic 
rays to magnetic fluctuations that have the 
same scale as the gyroradius of the c osmic 
ravs (iLerchei ll967tlK ulsrud & Pearcel 119691; 
I WentzelL 1 1 974I: ISkilUna 1 1 97 5 ) . When the per- 
turbed magnetic field grows to values similar 
to the background field, the resonance condi- 
tion is lost, and any magnetic field amplifica- 
tion beyond SB/B x 1 must be non-resonant. 
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With the development of diffusive shock 
acceleration, the need for efficient scattering 
around the shock region arised in order to get to 
cosmic ray energies of lO'^' eV. Additionally, 
observational evidence of narrow confinement 
regions around the shock for relativistic elec- 
trons required field strengths of the order of 
100 fjG, rather than just tens that could be 
expected based on compressi on of interstellar 
magnetic fields alone (e.g. IVink & Lamind, 
l2003HVolket~aLll2005h . The non-resonant in- 
stability developed by iBelll (|2004 resolved 
this. However, care should be taken by trans- 
lating a magnetic field strength inferred from 
observations of TeV electrons to the magnetic 
field strength confining the PeV cosmic rays. 

On scales smaller than the gyroradius ef- 
fective growh of magnetic fluctuations is gen- 
erated because of a plasma current. The cur- 
rent is a reaction to the cosmic-ray current, 
which requires balancing in the local plasma. 
Magnetic field fluctuations perpendicular to 
the current result in a J x B force that, for 
a given polarisation, focusses the cosmic rays 
within the growing magnetic loops, and pushes 
the plasma outwards. With a frozen-in mag- 
netic field, the stretching of the field lines 
results in growth. In the approximation that 
the larmor radius of the cosmic rays is much 
larger than the fluctuations, the cosmic ray 
current is effectively undisturbed and this can 
result in significant magnetic field amplifica- 
tion. Simulations of this instability in the non- 
linear regirne_Jndicate that this indeed is the 
case jBell 12004). An estimate for the satura- 
tion value of the magnetic field results in this 
value being proportional to the shock veloc- 
ity to the power 3/2. Therefore, especially for 
young remnants where the shock velocities are 
high, magnetic field growth to orders of mG 
seems to be feasible. 

The scales of the fluctuations triggered 
by this non-resonant instability are smaller 
than the gyro-radii of the driving cosmic rays. 
Scattering, being most efficient for cosmic rays 
having a larmor radius of the same order as 
the magnetic fluctuations, can become very 
efficient. Cosmic rays with low energies can 
therefore be confined near the shock region. 
However, in order to confine the high-energy 



end of the cosmic rays, magnetic fluctuations 
need to grow on a scale larger than the gyrora- 
dius of the dominant cosmic ray component. 

In a recent work, we describe the growth of 
magnetic fields on scales longer than the gyro- 
radius, due to the influence of the stress tensor 
arising fror n gradients in the pert urbed cosmic 
ray current (ISchure& Belli 1201 Ih . Small-scale 
turbulence, resulting from the non-resonant in- 
stability on small scales, induces an effective 
diffusivity that couples the right- and the left- 
hand polarisation through its action on the 
plasma. 

From th e disp ersion relation in 
ISchure & Bell (120111) . the influence of 
contributions from the stress tensor and col- 
lisions can be evaluated. The full dispersion 
relation, taking into account the stress tensor 
and collisions, looks like: 



2„2 



5(3v + icog) j 



(1) 



5(3v + icijg) ) 



The upper signs correspond to the left-handed 
polarisation (similar to the polarisation of a 
gyrating proton) and the lower signs to the 
right-hand polarisation (the one that is unsta- 
ble to the small-scale nonresonant Bell insta- 
bility). The magnetic fluctuations grow when 
the imaginary part of the square root of the 
above equation is positive, lj is the complex 
frequency, k the wavenumber, c speed of light, 
V the effective scattering frequency, a>g the gy- 
rofrequency, and Q. the growth rate of the non- 
resonant Bell instability that contains informa- 
tion on the driving current; D. - -sJkjoB/(pc), 
with B the unperturbed magnetic field strength 
parallel to k, jo the return current, and p the 
plasma mass density. 

The terms including factors of k^ result 
from the inclusion of the stress tensor Effects 
from the small-scale non-resonant instability 
are included through v, the effective scatter- 
ing frequency. Scattering on small scales is ex- 
pected to arise earlier than on long scales, since 
the growth rate increases rapidly for krg » 1 . 
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Bohm diffusion is the regime where the effec- 
tive scattering frequency is of the same order 
as the gyrofrequenc y, i.e. v ^ ajg. 

The method in ISchure & Bell ( 1201 lb re- 
covers the 2004 growth rate for the right po- 
larisation in the regime where kc » cOg. The 
resonant instability is not directly captured in 
the mono-energetic approach. When the col- 
lisionality is zero, v = 0, the only growing 
mode on small scales is the right-hand one, 
whereas on long scales only the left-hand mode 
is unstable. When v both modes are un- 
stable in the long-wavelength regime, although 
the left-hand mode continues to dominate until 
v/t Og > 1 / V3. 

iBvkov et al.l (1201 ih also investigated an in- 
stability on the large scales that takes into ac- 
count the effect of mean field dynamo the- 
ory. Averaging of the magnetic fluctuations on 
small scales, arising from the 2004 Bell in- 
stability, generates a ponderomotive force that 
acts to amplify the field on large scales. They 
found rapid growth, and their analysis also in- 
dicates that the left-hand mode is the dominant 
growing mode in the long-wavelength regime. 
A third option of growth on larger scales is that 
of growing small scale fluctuations. The fila- 
ments resulting from the small-scale nonreso- 
nant instabihty can m erge and grow in scales 
(iReville & Belli 1201 Ih . 



3. Typical growth times 

The driving current is a result of the stream- 
ing of cosmic rays, being proportional to the 
shock velocity: jo - n^^qv^, with «ci- the cos- 
mic ray number density, q the unit charge, 
and v, the shock velocity. Typical values for 
growth rates are summarized in Table [T] and 
plotted in Figure [T] A current density of jq - 
10"^ stA cm"^ is assumed at the shock, cor- 
responding to Her ~ 10"^ cm"^ and v, - 
lO** cm s-i. 

If Mcr is constant, the current scales with 

( ffv) ( 10') cm s-i )' gi^^" ^ powerlaw spectrum 
with a slope of p - 2, and a specific scat- 
tering efficiency. The particle energy is that 
of the driving cosmic rays, which is increas- 
ing with distance upstream of the shock when 
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Fig. 1. Growth time in years as a function of 
wavenumber. The solid (dashed) lines repre- 
sent the right- (left-) hand polarisation. Values 
are given for different efficiencies of the small- 
scale scattering rate. On long scales the left- 
hand mode has the quickest growth, whereas 
on scales smaller than resonance, the right- 
hand mode dominantes. 



diffusion propagates proportional to the par- 
ticle energy. Bohm diffusion can be consid- 
ered to be the minimal diffusion rate. Dividing 
this by the shock velocity gives a length- 
scale to which particles can diffuse upstream, 
given by L = D/v, = (rjc)/(3vj). The cor- 
responding time scale can be found by di- 
viding again over the shock velocity, to de- 
termine when the shock has caught up with 
the location where the instability is triggered, 
giving a lower limit of r = crg/(3v^) - 

0-35 (ifv) iM' iwik^y'^ y- The growth 
time scales with (^) (3^) {^j^) , 
for Her - 10"^. The number of e-foldings for 
growth, jT, can be seen to scale with Vj^^^, 
leaving only the shock velocity as a free pa- 
rameter to increase the linear growth effec- 
tively. 

Alternatively, it can be considered that ric,-, 
and therefore jo, scales with the kinetic en- 
ergy of the shock, such that a fixed propor- 
tion of pvj is converted into cosmic rays. For 
a mono-energetic approach, this translates to 
TicrEcr - Tjjpvl such that the current den- 
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Table 1. Linear growth times for the right (and left between brackets) hand modes, given 
for a particle energy of 1 TeV, magnetic field strength B - 3 fiG, shock velocity - 
lO^cm s"\ and current jo - 10"^ stA cm"^. It can be derived for other values by us- 
ing the fact that it scales with {^^^iy-^^ ( lo** cm s-' ) constant in time, or with 
(rfv)(3^) ' ( wo cm s-' ) ^^^(oom) '^^ if "cr is a fraction of the shock kinetic energy. 





r ' (yr) 


7 ' (yr) 


r ' (yr) 


7 ' (yr) 


J-' (yr) 


7 ' (yr) 




(fcrj = 0.1) 


{krg = 0.5) 


{krg = 0.8) 


(kr, = 1) 


(kr, = 5) 


(krg = 10) 


1/10 


600 (119) 


50.2(10.5) 


22.1 (5.00) 


4.47 (1.10) 


0.651 (13.1) 


0.498 (47.4) 


1/3 


264 (164) 


23.0(14.7) 


10.8 (7.31) 


2.34(1.67) 


0.700 (4.62) 


0.503 (14.6) 


1/V3 


245 (246) 


21.7 (22.2) 


10.6 (11.2) 


2.35 (2.58) 


0.776 (3.53) 


0.513 (8.92) 


1 


283 (462) 


25.3 (41.6) 


12.5 (20.9) 


8.95(15.2) 


0.949 (3.58) 


0.542 (6.01) 



sity becomes jo - ricrqvs - riq^pv^jEcr 
and thus scales with (^) ' (1057^^)^ The 
growth time in that case scales as y"' oc 

( tI?) (3^)"' (to^)"'^' {muV^ and the 
number of e-foldings is proportional to 

— 1/2 

yr oc V5 . Of course in these equations 

{^){^) ' can be replaced by ( 
with Kg the gyroradius. The numbers in Table[T] 
correspond to a value of 77 = 0.014. It is ex- 
pected that the fraction of kinetic energy go- 
ing into cosmic rays can be significantly higher 
than this, resulting in shorter growth times by 
a factor of the square root of this efficiency. 
For an efficiency of 10% this would result in 
reduced growth times compared to those in 
Table[T]by a factor 2.7. 

4. Discussion and conclusion 

Acceleration of cosmic rays in supernova rem- 
nants requires the growth of magnetic fluctu- 
ations in order to confine cosmic rays to the 
shock region. For the highest-energy cosmic 
rays the non-resonant regime on scales larger 
than the gyroradius of the dominant cosmic ray 
component is of importance in confining cos- 
mic rays as they reach higher energies. A non- 
resonant instability is discussed that is driven 
by the return current in the plasma, result- 
ing from a reaction to ffie cosmic ray cur- 



rent. The stress-tensor is crucial in generating 
an instability in the long-wavelength regime. 
Scattering on small scales aids in mediating the 
instability between the two polarisations. The 
number of e-folding times for the instability at 
a given scattering efficiency is shown to depend 
only on the shock velocity, and typical growth 
times are given. 
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